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THEORETICAL  AND  EXPERIMENTAL  INVESTIGATION  OF  THE 

EFFECT  OF  TUNNEL  WALLS  ON  THE  FORCES  ON  AN 

OSCILLATING  AIRFOIL  IN  TWO-DIMENSIONAL 

SUBSONIC  COMPRESSIBLE  FLOW1 

By  Harry  L.  Runyan,  Donald  S.  Woolston, 
and  A.  Gerald  Rainey 

SUMMARY 


This  paper  presents  a  theoretical  and  experimental  investigation 
of  the  effect  of  wind-tunnel  walls  on  the  air  forces  on  an  oscillating 
wing  in  two-dimensional  subsonic  compressible  flow.  A  method  of  solving 
an  integral  equation  which  relates  the  downwash  on  a  wing  to  the  unknown 
loading  is  given,  and  some  comparisons  are  made  between  the  theoretical 
results  and  the  experimental  results.  A  resonance  condition,  which  was 
predicted  by  theory  in  a  previous  report  (NACA  Rep.  1150),  is  shown 
experimentally  to  exist.  In  addition,  application  of  the  analysis  is 
made  to  a  number  of  examples  in  order  to  illustrate  the  influence  of 
walls  due  to  variations  in  frequency  of  oscillation,  Mach  number,  and 
ratio  of  tunnel  height  to  wing  semichord. 


INTRODUCTION 


In  the  evaluation  of  results  obtained  by  measurement  of  the  forces 
on  a  wing  in  a  wind  tunnel,  the  question  of  the  effect  of  the  tunnel 
walls  arises.  In  the  case  of  steady  flow  the  problem  has  been  extensively 
investigated  and,  in  general,  relatively  simple  factors  have  been  deter¬ 
mined  which  can  be  used  to  modify  measurements  of  the  forces  on  a  wing 
in  a  tunnel  to  correspond  to  free-air  conditions.  However,  the  corre¬ 
sponding  problem  of  the  effect  of  walls  on  an  oscillating  airfoil  has 
received  relatively  little  attention,  particularly  in  the  case  of  com¬ 
pressible  flow.  The  present  paper  concerns  the  wall  effects  in  the 
oscillating  case  and  treats  the  problem  in  two-dimensional  subsonic 
compressible  flow. 

In  incompressible  flow,  theoretical  treatments  of  wall  effects  on 
oscillating  wings  have  been  made  by  several  investigators  and  reported 
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in  references  1,  2,  and  3.  These  investigators  have  shown  generally 
that  the  tunnel- wall  effects  are  a  maximum  for  some  small  values  of  the 
reduced  frequency  and  that  the  wall  effects  become  negligible  as  the 
reduced  frequency  is  increased.  Extension  of  the  theoretical  treatment 
of  the  problem  to  include  the  effects  of  compressibility  of  the  fluid 
has  been  reported  in  reference  4.  In  this  reference,  it  is  shown  that, 
in  addition  to  the  large  effect  noted  at  low  values  of  the  reduced  fre¬ 
quency,  under  certain  conditions,  large  effects  of  the  walls  may  be 
encountered  at  higher  values  of  the  reduced  frequency.  These  effects 
are  due  to  an  acoustic  resonant  phenomenon  which  occurs  when  a  disturbance 
from  the  oscillating  wing  is  reflected  from  the  tunnel  wall  back  to  the 
wing  with  such  a  phase  relationship  that  it  reinforces  a  succeeding 
disturbance . 

In  reference  4  the  problem  was  expressed  as  an  integral  equation 
which  relates  the  known  downwash  distribution  over  the  airfoil  to  the 
unknown  lift  distribution.  One  purpose  of  the  present  paper  is  to  dis¬ 
cuss  further  the  integral  equation  and  to  demonstrate  a  method  of  solving 
it.  A  second  purpose  is  to  present  some  results  showing  wall  effects 
calculated  by  this  procedure  and,  in  some  cases,  to  compare  the  calcu¬ 
lated  results  with  experimental  results.  This  phase  of  the  investigation 
is  given  in  three  parts:  (l)  a  comparison  between  analytically  and  exper¬ 
imentally  determined  values  for  the  lift  and  moment  on  a  wing  oscillating 
in  pitch  at  several  subsonic  Mach  numbers;  (2)  an  analytical  study  of  the 
effects  of  a  variation  in  Mach  number  for  a  constant  ratio  of  tunnel  height 
to  wing  semichord;  and  (3)  an  analytical  study  of  the  effects  of  a  vari¬ 
ation  in  the  ratio  of  tunnel  height  to  wing  semi chord.  Portions  of  this 
material  have  been  reported  previously  in  reference  5  and  are  included  in 
the  present  paper  in  order  to  provide  a  more  extensive  and  unified 
presentation. 

As  a  check,  the  integral  equation  for  the  downwash  on  a  wing  oscil¬ 
lating  between  walls  in  a  compressible  medium  is  reduced  to  the  zero- 
frequency  condition  and  is  shown  in  an  appendix.  The  resulting  expression 
is  in  agreement  with  steady-state  results. 

The  calculation  procedure  and  the  results  contained  in  this  paper  are 
of  significance  for  such  problems  as  the  experimental  measurement  of  the 
forces  on  an  oscillating  airfoil,  the  determination  of  wing-flutter  char¬ 
acteristics  in  wind  tunnels,  and  also  in  certain  possible  types  of  flutter 
of  airfoils  in  cascade. 


SYMBOLS 


a  velocity  of  sound,  ft/sec 

An  coefficients  in  series  expression  for  lift  distri¬ 

bution  (eq.  (l6))  ,  where  n  =  0,  1,  2,  .  .  • 
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0  =  -cos^^x 
eQ  =  -cos-^Xq 


fluid  density,  slugs/cu  ft 

phase  angle  between  lift  force  and  position  of 
pitching  wing,  deg 

phase  angle  between  lift  force  and  position  of 
translating  wing,  deg 

phase  angle  between  moment  and  position  of  pitching 
wing,  deg 

phase  angle  between  moment  and  position  of  trans¬ 
lating  wing,  deg 

circular  frequency  of  oscillation,  radians/sec 

circular  frequency  at  resonance,  radians/sec 

pressure  difference  between  upper  and  lower  surface, 
lb/sq  ft 

Primed  quantities  refer  to  a  wing  in  free  air. 
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This  section  is  concerned  with  the  development  of  a  method  for 
solving  the  integral  equation,  originally  derived  in  reference  4,  which 
relates  the  downwash  to  the  loading  on  an  oscillating  wing.  The  basic 
integral  equation  and  its  kernel  is  given  by  equations  (l)  and  (2) . 
Reduction  of  the  kernel  is  made  in  equations  (3)  to  (10) .  Alternative 
series  expressions  for  the  kernel  which  are  suitable  for  numerical 
computation  are  given  by  equations  (ll)  to  (15)  •  The  loading  on  the 
wing  is  given  by  equation  (l6),  the  downwash  expression  by  equations  (l8) 
and  (19)>  and  finally  the  lift  and  moment  expressions  by  equations  (20). 
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The  Integral  Equation  and  Its  Kernel  Function 

The  integral  equation. -  The  integral  equation  of  reference  4  for 
the  vertical  velocity  or  downwash  of  an  oscillating  airfoil  between 
plane  walls  may  be  written  as 


w(x) 


L(xq)[k(M,z)  +  K(M,z,H)]  dxQ 


(1) 


where  w(x)  is  the  known  vertical  velocity  (or  known  motion  of  the  wing) 
and  L(x0)  is  the  unknown  lift  distribution  or  the  local  strength  of  a 
distribution  of  oscillating  pressure  doublets.  The  functions  within  the 
brackets  comprise  the  kernel  function  of  the  integral  equation  and  appear 
formally  as 


K(M, z)  =  T-i-  lim  e-ik(x-x0) 
y — »0 


'x“xo  T&S 


°eP2  &  +  (2a) 


ik 

K(M,Z,H)  =  rd—  lim  e"ik(x-xo)  P  °  2  (-l)n  eP2?  |!L  He/2) 


-oo  n=l 


®  \ft 2  +  32(y  -  nH)2 
P2 


d|  (2b) 


The  first  function  K(M,  z)  corresponds  to  the  kernel  for  the  free-air 
condition  as  given  by  Possio  (ref.  6).  The  second  function  K(M,z,H), 
containing  the  infinite  summation,  is  the  additional  part  of  the  kernel 
arising  from  the  effect  of  the  walls.  Physically,  a  kernel  function 
represents  the  contribution  to  the  vertical  velocity  at  a  field  point 
due  to  a  pulsating  pressure  doublet  of  unit  strength  located  at  any  other 
point  in  the  field.  For  the  particular  case  represented  by  equations  (2), 
the  kernel  function  gives  the  vertical  velocity  in  the  plane  of  a  wing 
located  in  the  center  of  the  tunnel.  The  expression  K(M,z)  gives  the 
downwash  of  a  doublet  in  the  plane  of  the  wing,  whereas  the  expression 
K(M,z,H)  gives  the  downwash  due  to  the  system  of  images  which  mathemat¬ 
ically  represents  the  walls. 

Reduction  of  the  kernel  function. -  The  integrals  contained  in  the 
expressions  for  the  kernel  function  in  equations  (2)  are  improper  because 
they  have  an  infinite  limit  and  also  because,  at  certain  points,  the 
integrands  become  singular.  This  section  is  concerned  with  the  reduction 
of  these  integrals  to  a  form  more  amenable  to  computation. 

By  making  use  of  the  fact  that  the  Hankel  functions  in  equations  (2) 
satisfy  the  identity 
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The  integrals  of  equation  (1)-)  that  contain  partial  derivatives  of 
Hankel  functions  can  he  integrated  twice  by  parts  to  obtain 


“°  ^  ft  HO* 

d!2 


(2)  Mk 
P2 


£2  +  .P2(y  -  nH)2  <*t  = 


* 
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The  last  integral  of  equation  (5)  may  he  written  in  two  parts  as 


The  first  integral  on  the  right-hand  side  of  equation  (6)  will  be  left 
temporarily  in  integral  form  and  will  he  treated  in  the  following  section. 
(See  evaluation  of  S3  following  eq.  (13).) 

The  second  integral  on  the  right-hand  side  of  equation  (6)  has  not 
been  integrated  in  closed  form;  however,  in  wind-tunnel  problems  it  can 
be  handled  conveniently  by  approximate  methods.  (An  alternative  means 
of  treating  this  integral,  which  avoids  the  approximation,  but  is  somewhat 
more  tedious,  will  be  indicated  in  the  section  following  eq.  14(c).) 

A  practical  assumption  which  is  often  made  in  the  analysis  of  the  effect 
of  wind-tunnel  walls  is  that  the  tunnel  height  is  considered  large  compared 
with  the  wing  semichord.  With  this  assumption  the  argument  of  the  Hankel 
function  in  equation  (6)  can  be  written  as  (in  the  limit  as  y— *0) 


provided  that  — «  1. 

pnH 

This  approximation  implies  that  the  airfoil  images,  and,  in  partic¬ 
ular,  the  closest  image  n  =  1,  are  a  sufficient  distance  from  the  airfoil 

so  that  the  actual  distance  \f^  +  p^(nH)^  may  be  replaced  by  the  vertical 
distance  |3nH  of  the  image  above  the  airfoil.  Of  course,  this  approxima¬ 
tion  does  not  hold  for  Mach  numbers  close  to  or  equal  to  unity.  With  this 
approximation,  the  second  integral  of  equation  (6)  can  be  expressed  as 
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lim 
y— >0  ~o 


-xD  00  . 

eP2  (-l)”  %(2) 

n=l 


*§  /l2  +  P2(y  -  nH)2 
P2 


-?  V 

MknM  /'x"xo  r#l 


(-DnHo (2)(*=)  /  %P2Sd| 

n=l  V  P  /  J0 


■s:(-«n^(2)(T)^ 


[\£: 

J(MB)  -  1 

/ik 

(7) 


and  these  equations  may  be  used  to  express  equation  (4)  as 

■1 


v(x)  = 


oib 

pU2 


— 

L(x0)  K(M,z) 
1  L 


+  K(M,z,H) 


dx, 


(8) 


where 


K(M,z)  =  e~iz  1  eP2 


iz  r- 


-Hq^HpRo)  +  ■1^~— ~  Hi^2)(pRo) 


IP' 


—  loge  —  +  -  +  /  M  eiu  Ho^2)  (M|u|)du 
tc  p  M  J  q 


and 


K(M,z,H)  =  ^ 


-iz 


1Z 


P2  V  (-l)n  Ho(2)(pRn)  + 


n=l 


(ep2  -  l)  'll  (-Dn  H0<2)(M) 


(9) 


>00  ik, 


ik  >  (-l)n  /  e  P2^  Hq^2) 

n=l  J0 


M  y/|2  +  p2(nH)2 

P2 


d|  + 


e?  f:  (.1)”  1M(X  -  *°.>  Hl(2) 


n=l 


(pRn) 


(10) 
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in  which  use  has  been  made  of 

n  =  u  =  v  = 

p2  p25  p2 

R0  =  |x  -  xQ|  Rn  =  \[(x  -  xQ)2  +  p2(nH)2 

Equation  (8),  together  with  the  definition  of  equations  (9)  and  (10), 
permits  the  determination  of  the  effect  of  tunnel  walls  on  a  lift  distri¬ 
bution  L(0O)  for  a  given  downwash  distribution  w(x) .  The  integral 
equation  for  the  case  of  no  tunnel  walls  checks  the  results  of  Possio 
(ref.  6).  For  the  case  with  walls  and  for  the  limiting  steady-flow  case 
of  zero  frequency,  it  is  possible  to  obtain  a  mathematical  check  with 
some  existing  results;  this  check  is  shown  in  the  appendix. 

Alternative  series  expressions  for  kernel.-  Although  the  form  of  the 
kernel  K(M,z,H),  given  by  equation  (10),  could  be  used  for  calculation, 
alternative  series  which  are  more  highly  convergent  may  be  used  and  are 
given  in  this  section. 

The  kernel  K(M,  z,H)  is  the  sum  of  four  infinite  series  which  can 
be  written  as 

— i  z 

K(M,z,H)  =  +  C2S2  +  C5S5  +  (11) 

where  the  Sn's  denote  the  indicated  infinite  summations  of  equation  (10) 
and  the  Cn*s  the  respective  multipliers. 

Series  Sp  and  S2  of  equation  (ll)  may  be  put  in  a  more  rapidly 

convergent  form  according  to  Infeld,  Smith,  and  Chien  (ref.  7)  •  When  the 
variables  p  and  e  are  introduced,  where 

p  =  MkH 
2*13 

and 

x  -  x0 

the  series  S^  and  S2  can  be  written  as 
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and. 


L*J 

Si  =  21  (-1)n  Ho(2)(nRn) 
n=l 


n=l 


(~l)n  Hq(^^2jip  l/e2  +  n2  j 


2ie 


-rte  \|l-4p2 


_  _1 _ 

2jt^  v^1  -  4p2 


+  21 


e  J 


/(2n+l)2-Ap2 


n=l  ^/(2n  +  1)2  -  4p2 


+ 


.  ,Ho<2)(2Mp) 

\/(2n  -  l)2  -  4p2 


(12) 


S2  =  2Z  (-1)11  Hq(2) 


(-l)n  Hq(  2)  ( 2rtnp) 


n=l 


_1_ 

2it 


00 

-it  +  21(  y  +  log~  +  hi.  / 

\ 

1  1  1 

n=l 

2n  -  1)^  -  2n  -  1 

(13) 


where  Euler's  constant  7  =  0.577215* 

Series  S5  may  "be  evaluated  "by  utilizing  the  expression  for  Si 
(eq.  (12))  and  integrating  the  resulting  expression  to  obtain 


s3  = 


ikt 

r  .-p* 

00 

r  “i 

(-l)nHo<2) 

+  p2(nH)  2 

1 0 

n=l 

_p 

d| 


(14a) 


y 
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d| 


(14b) 


S5  =  -  J-  logg  ^£-8.  + 
■>  rtk  ***  M 


2n  +  1)2  _  4P2 

(14c) 

It  is  of  interest  to  note  that  series  Sj  may  he  employed  in  an 

alternative  means  of  integrating  equation  (7).  For  application  to  wind- 
tunnel  problems  ,  where  the  ratio  of  tunnel  height  to  wing  semichord  is 
small,  or  in  application  to  cascade  problems,  the  approximation  employed 
in  integrating  equation  (7)  becomes  less  valid.  It  is  possible  to  avoid 
the  use  of  the  approximation  by  writing  the  integral  of  equation  (7)  in 
a  form  which  is  identical  to  that  of  equations  (l4a)  and  (iVb)  with  the 
exception  of  the  upper  limit.  The  integral  containing  the  Hankel  function 
can  be  evaluated  by  employing  the  tables  of  Schwarz  (ref.  8).  The  second 
integral,  containing  only  an  exponential  term,  can  be  integrated  in 
closed  form,  as  was  done  to  obtain  equation  (l4c)  . 


Series  Sl^_  may  be  evaluated  in  a  direct  manner  by  employing  tables 
of  the  Hankel  function  and  by  using  for  large  values  of  the  argument  the 
appr  oximat i on 

With  the  aid  of  series  Sj_,  Sq,  Sj,  and  S4,  the  kernel  K(M, z,H)  may 
be  evaluated. 
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Method  of  Solution 

A  method  of  using  equation  (8)  to  determine  the  aerodynamic  forces 
on  a  wing  oscillating  in  the  presence  of  plane  walls  is  now  discussed. 

The  method  under  consideration  is  one  of  collocation  similar  to  that 
used  by  Possio  (ref.  6)  and  Frazer  (ref.  9)  for  the  case  of  no  walls. 

The  approach  involves  the  assumption  of  an  appropriate  series  expression 
for  the  lift  distribution,  substitution  of  this  series  in  the  integral 
equation  for  the  downwash,  and  calculation  of  the  downwash  at  arbitrarily 
selected  points  on  the  chord  (control  points).  Thus  equation  (8)  is 
reduced  to  a  set  of  simultaneous  equations,  the  unknowns  of  which  are 
the  coefficients  of  the  assumed  expression  for  the  loading. 

Expression  for  the  loading.-  The  expression  which  is  assumed  for  the 
lift  distribution  is  a  trigonometric  series  expansion  which  satisfies  the 
Kutta  condition  at  the  trailing  edge  and  which  has  the  proper  type  of 
singularity  at  the  leading  edge.  This  expression  is 


00 

HiEgl  =  Aq  cot  ^  An  sin  n0o  =  L(90)  (l6) 

plK  ^  n=l 

where  xQ  =  -cos  0O  and  the  AnTs  are  unknown  coefficients  to  be  deter¬ 
mined  in  accordance  with  the  downwash  w(x) ,  which  is  known  from  the 
motion  of  the  wing.  It  is  desirable  to  rewrite  equation  (8)  in  terms  of 
the  variable  0O  as  follows: 


w(x)  =  Uk 


L(0O)  K(M,  z)  sin  0O  d0o  + 


L(0O)  K(M,z,H)  sin  0O  d0o 


(IT) 


The  first  integral  on  the  right-hand  side  of  equation  (17)  is  the 
integral  expression  first  derived  by  Possio  (ref.  6)  for  the  condition 
of  no  walls.  Its  solution  has  been  treated  by  several  investigators  and 
will  not  be  discussed  here.  It  can  be  expressed  entirely  in  terms  of 
the  unknown  coefficients  An  of  equation  (l6) .  The  second  integral  of 
equation  (17)  may  be  evaluated  by  the  use  of  equations  (12),  (13),  (l4), 
and  ( 15) . 


Determination  of  the  aerodynamic  forces. -  The  integrals  of  equa¬ 
tion  (17) are  determined  for  a  selected  number  of  control  points  and 
equated  to  the  expression  for  the  downwash.  The  expression  relating  the 
downwash  to  the  motion  of  a  wing  translating  (h)  and  pitching  (a)  about 
an  axis  located  at  xa  is 


* 


w(x)  =  h  +  Ua  +  b(x  -  xa)a 


(18) 
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or,  with  the  assumption  of  harmonic  motion. 


u 


*a) 


a 


(19) 


Equation  (19)  is  used  to  calculate  w(x)  for  values  of  x  appropriate 
to  each  of  the  selected  control  points.  A  set  of  simultaneous  equations 
can  then  he  written,  the  number  of  which  corresponds  to  the  number  of 
control  points  employed  and  (conveniently)  to  the  number  of  terms 
retained  in  the  series  for  L(0O).  The  unknown  coefficients  may  now  be 
determined  by  solving  these  simultaneous  equations.  The  total  lift  and 
moment  about  the  midchord  are  given  in  terms  of  the  coefficients  An 
through  the  relations 


Effect  of  the  number  of  control  points  considered.-  An  investigation 
was  made  of  the  number  of  terms  of  the  series  for  the  lift  distribution 
(eq.  (l6))  and  thus  of  the  number  of  control  points  required  to  obtain 
satisfactory  accuracy.  Calculations  were  performed  for  a  particular  case 
by  increasing  the  number  of  control  points  and  the  number  of  terms  of  the 
loading  series  until  the  solutions  were  in  reasonable  agreement.  For  the 
case  considered,  three  terms  of  the  series  for  the  lift  and  three  control 
points  at  the  quarter-,  half-,  and  three -quarter -chord  positions  gave 
satisfactory  results.  The  consideration  of  two  additional  control  points 
at  the  leading  and  trailing  edges,  together  with  two  additional  terms  of 
the  lift  series,  made  no  significant  change  in  the  results.  For  high 
values  of  the  reduced-frequency  parameter  k,  the  use  of  additional 
control  points  might  be  necessary. 


The  procedure  just  discussed  involves  consideration  of  a  continuous 
distribution  of  pressure  doublets  over  the  chord.  Calculations  requiring 
much  less  computing  can  be  made  by  considering  the  chordwise  loading  to 
be  concentrated  in  a  single  doublet  located  at  the  quarter  chord  and 
satisfying  the  downwash  at  the  three-quarter  chord.  In  the  case  of  the 
lift,  this  approach  has  been  found  to  give  fairly  good  agreement  with 
the  results  of  the  more  elaborate  calculations  except  in  the  vicinity 
of  the  resonant  frequency. 
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The  Analytically  Indicated  Resonance  Phenomenon 

Two-dimensional  tunnel. -  By  examination  of  equations  (12)  and  (13),  , 

it  may  be  seen  that  the  series  Sq  and  S2  become  infinite  when 

hv2  =  (2n  -  l)2 

or  where 

^  =  jt3(2n  -  1)  (n  =  1,2,3,  •  •  •)  (2l) 

a 

At  these  critical  values  of  the  frequency  parameter,  the  expression  for 
the  kernel  K(M,z,H)  (eq.  (ll))  becomes  infinite  for  all  values  of  x. 
Physically,  this  condition  represents  a  resonance  in  the  tunnel  involving 
a  transverse  oscillation  of  the  moving  air  between  the  walls. 

The  fundamental  or  smallest  critical  values  of  o>H/a  corresponding 
to  n  =  1  in  equation  (21)  are  shown  plotted  as  functions  of  Mach 
number  M  in  figure  1.  Equation  (21)  and  figure  1  show  that  finite  0 

values  of  the  critical  frequency  exist  for  the  condition  M  =  0,  U  =  0, 
and  a  ^  oo.  These  conditions  correspond  to  a  compressible  fluid  at  zero 
velocity  in  the  tunnel.  As  the  Mach  number  is  increased,  the  critical-  * 

frequency  parameter  decreases  rapidly  and  becomes  zero  at  a  Mach  number 
of  unity. 

As  indicated  by  equation  (l),  the  product  of  the  lift  and  the  kernel 
function  must  remain  equal  to  the  vertical  velocity  over  the  wing}  this 
velocity  is  defined  by  the  motion  of  the  wing  and  remains  finite.  The 
product  of  the  lift  and  the  kernel  function  can  remain  finite  only  if  the 
lift  approaches  zero  as  the  kernel  becomes  infinite.  This  condition  in 
the  tunnel  is  analogous  to  the  well-known  case  of  a  simple  undamped- 
spring-mass  system  for  which,  at  the  resonant  frequency,  theory  predicts 
an  infinite  deflection  of  the  mass  occurring  even  with  a  forcing  function 
of  small  amplitude. 

Circular  tunnel.-  A  resonance  can  also  be  demonstrated  for  the 
infinite  circular  tunnel.  The  nature  of  the  boundary- value  problem,  for 
this  case,  makes  it  possible  to  separate  variables;  therefore,  the 
governing  partial-differential  equation  can  be  reduced  to  Bessel's 
equation.  (See,  for  instance,  ref.  10.)  The  resonant  frequencies  are 
then  found  as  the  roots  of  the  equation 

Jn'(i?p)  =  0  (n  =  0'1,2,  *  ' 

w 
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or 


mD 

2a 


PnP 


where  Jn  represents  the  Bessel  function  of  the  first  kind,  D  is 
the  tunnel  diameter,  and  is  the  root  of  the  equation 

Jn'(Pn)  =  0 


Values  for  pjj  for  the  first  several  modes  are  =  1.84,  3.05,  and 
4.17-  Note  that,  for  a  circular  tunnel  having  a  diameter  equal  to  the 
height  of  a  plane  tunnel,  the  fundamental  frequency  is  3.68/rt  =  1.17 
higher  than  resonant  frequency  in  the  plane  tunnel  discussed  in  this 
paper . 


EXPERIMENTAL  INVESTIGATION 
Wind  Tunnel 


The  experimental  part  of  the  investigation  of  the  effect  of  tunnel 
walls  on  the  forces  acting  on  an  oscillating  airfoil  was  conducted  in 
the  Langley  2-  by  4-foot  flutter  research  tunnel.  For  these  tests,  a 
rectangular  test  section  having  dimensions  of  2  feet  by  3.8  feet  was 
used.  This  tunnel  is  of  the  closed-throat,  single-return  type  and 
employs  either  air  or  Freon- 12  as  a  testing  medium  at  pressures  from 
1  atmosphere  down  to  about  1/8  atmosphere. 

It  has  been  shown  previously  that  the  resonant  frequency  varies 
directly  as  the  speed  of  sound.  Inasmuch  as  Freon-12  has  a  speed  of 
sound  equal  to  about  one-half  that  of  air,  the  experiments  to  be  dis¬ 
cussed  were  conducted  in  Freon-12  so  that  the  resonant  frequency  could 
be  surveyed  within  the  frequency  limitations  of  the  equipment. 


Model  and  Oscillating  Mechanism 

Figure  2  is  a  schematic  drawing  of  the  test  section  with  the  model 
and  oscillating  mechanism  installed.  The  model  had  a  chord  of  1  foot 
and  MCA  65-OIO  airfoil  sections)  it  completely  spanned  the  2-foot  dimen 
sion  of  the  test  section.  The  gaps  between  the  model  and  the  tunnel 
wall  were  sealed  by  end  plates  which  rotated  with  the  model.  The  model, 
driven  symmetrically  from  both  ends,  was  oscillated  in  pitch  about  the 
midchord  by  a  direct-drive  eccentric-cam  system  powered  by  an  induction 
motor  with  variable  frequency  supply. 


Instrumentation 


The  lift  and  moment  on  the  wing  were  obtained  by  electrical  integra¬ 
tion  of  the  outputs  of  12  model  49-TP  NACA  miniature  electrical  pressure 
gages.  The  pressure  gages,  which  are  described  in  considerable  detail 
in  reference  11,  were  located  at  the  center  of  the  span  at  2-5>  10, 

15,  20,  50,  40,  50,  60,  70,  80,  and  90  percent  of  the  chord.  Each  gage 
was  arranged  to  indicate  the  difference  in  pressure  between  orifices  on 
the  upper  and  lower  surfaces.  Electrical  integration  techniques  used  in 
these  experiments  are  discussed  in  reference  12.  The  so-called  square- 
wave  method  of  weighting  was  used)  that  is,  the  pressure  indicated  by 
each  gage  was  assumed  to  represent  the  pressure  acting  over  a  portion  of 
the  chord  extending  one-half  the  distance  to  the  next  gage  both  forward 
and  rearward.  For  example,  the  fraction  of  the  chord  assigned  to  the 
first  gage  was  3-75  percent  and  to  the  sixth  gage  was  10  percent.  Some 
of  the  implications  of  this  method  of  integration  will  be  discussed  in 
a  subsequent  section. 

The  angular  displacement  at  the  midspan  position  was  indicated  by 
resistance-wire  strain  gages  attached  to  a  torque  rod  running  through 
the  center  of  the  hollow  wing.  One  end  of  the  torque  rod  was  fixed  to 
the  center  of  the  wing  and  the  other  end  was  fixed  to  the  tunnel  wall. 

A  schematic  diagram  of  the  instrumentation  is  shown  in  figure  3* 

The  magnitude  of  the  vector  representing  the  fundamental  component  of 
lift  or  moment  and  angular  displacement  was  indicated  on  an  alternating- 
current  vacuum-tube  voltmeter  attached  to  the  output  of  a  variable- 
frequency,  narrow-pass -band  filter.  In  essence,  the  filter  performed 
the  function  of  a  Fourier  analysis  in  that  both  random  components  and 
higher  harmonics  were  removed  from  the  signal.  In  order  to  measure  the 
phase  angle  between  lift  or  moment  and  the  angular  displacement,  the 
output  of  the  filter  was  fed  into  a  pulse-shaping  circuit  designed  to 
convert  the  sinusoidal  signals  into  pulses  corresponding  in  time  to  the 
"crossover"  points  of  the  original  signal.  The  pulses  were  then  fed 
into  an  electronic  chronograph  that  accurately  indicated  the  time 
interval  between  the  leading  pulse  which  started  the  chronograph  and 
the  lagging  pulse  which  stopped  it.  The  ratio  of  this  time  interval 
to  the  period  of  oscillation,  when  multiplied  by  360°,  yields  the  phase 
angle  in  degrees.  The  period  and  frequency  of  oscillation  were  deter¬ 
mined  by  starting  and  stopping  the  chronograph  with  the  angular- 
displacement  signal.  In  order  to  minimize  the  effects  of  small  differ¬ 
ences  in  components  between  the  two  circuits,  a  "tare"  switch  was  pro¬ 
vided  which  fed  a  single  signal  (the  angular  displacement)  through  both 
circuits.  The  resulting  phase  angle  represented  the  phase  shift  intro¬ 
duced  by  the  filters  and  pulse-shaping  circuits. 
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Test  Conditions 

The  Mach  number  of  the  tests  was  varied  from  M  =  0.35  to  M  =  0.70 
and  the  Reynolds  number  was  held  constant  at  about  5  X  10°  by  varying  the 
density.  The  frequency  of  oscillation  was  varied  from  0  to  60  cycles  per 
second,  and  the  magnitude  of  angular  displacement  was  about  1.2°  except 
for  some  lift  data  at  M  =  0.71  which  was  obtained  at  an  angular  dis¬ 
placement  of  about  2.4°. 


DISCUSSION  OF  RESULTS 


The  theory  and  calculation  procedure  and  the  experimental  technique 
discussed  previously  for  the  determination  of  the  forces  acting  on  a 
wing  oscillating  between  walls  have  been  applied  to  a  number  of  specific 
examples.  The  investigation  has  been  divided  into  three  parts:  (l)  A 
comparison  is  made  Of  analytical  and  experimental  results  obtained  for 
the  lift  and  moment  on  a  pitching  wing  for  several  subsonic  Mach  numbers, 
(2)  theoretical  results  for  the  effects  of  a  variation  in  Mach  number  at 
constant  tunnel  height  are  given  for  a  pitching  wing  and  also  for  a  wing 
undergoing  vertical  translation,  and  (3)  theoretical  results  for  the 
effects  of  a  variation  in  the  ratio  of  tunnel  height  to  wing  semichord 
are  presented  for  particular  values  of  Mach  number. 


Comparison  of  Theory  and  Experiment 

In  figure  4  a  comparison  is  made  of  analytical  and  experimental 
results  for  a  wing  oscillating  in  pitch  about  its  midchord.  Figures  4(a), 
4(b),  4(c),  and  4(d)  apply,  respectively,  to  Mach  numbers  of  0.35>  0.5, 
0.6,  and  0.7.  The  results  apply  to  a  ratio  of  tunnel  height  to  wing 
semichord  H  of  7-60. 

The  plots  on  the  left-hand  side  of  each  figure  show  the  magnitudes 
of  the  forces  and  moments  as  a  function  of  the  frequency  of  oscillation 
whereas  those  on  the  right-hand  side  show  the  corresponding  phase  angles. 
The  magnitudes  are  presented  as  ordinates  in  the  form  of  ratios  |  I^/l^  *  | 
and  |  Hx 1 1  •  In  these  ratios,  the  quantities  1^  and  are,  respec¬ 

tively,  the  lift  force  and  the  moment  on  a  wing  in  a  tunnel;  1^'  and 
1^'  are  the  theoretical  lift  and  the  theoretical  moment  on  a  wing  in 
free  air.  The  effect  of  the  tunnel  walls  appears,  therefore,  as  a  devia¬ 
tion  from  unity  of  the  ratios  |  |  and  |  when  1^  and 

Hx  are  the  theoretically  derived  forces  and  moments.  When  Iq,  and 
represent  the  experimental  forces  and  moments,  the  deviation  from  unity 
may  not  be  completely  attributed  to  the  effect  of  tunnel  walls  because 
such  factors  as  airfoil  thickness  and  viscosity  may  cause  deviation  from 
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the  elementary  theory.  The  abscissa  in  the  figures  is  the  ratio  of  the 
frequency  of  the  pitching  oscillation  to  a  frequency  calculated  for  the 
resonant  condition. 

Excellent  agreement  between  theory  and  experiment  is  obtained  for 
the  phase  angles,  in  most  cases,  for  both  the  lift  and  the  moment. 
Quantitatively,  however,  the  agreement  between  theory  and  experiment  for 
the  magnitudes  of  the  forces  is  not  as  good,  although  very  similar  trends 
are  demonstrated;  in  most  cases,  a  systematic  difference  appears.  Some 
possible  sources  of  the  differences  between  theory  and  experiment  are 
discussed  in  the  following  section. 

Examination  of  figure  4  reveals  that  the  theory  predicted  the 
resonant  frequency  very  well.  In  all  cases,  the  minimum  lift  and 
moment  was  found  to  lie  very  close  to  the  analytically  indicated  resonant 
frequency.  Theoretically,  the  lift  and  moment  reduce  to  zero  at  the 
resonant  condition.  Under  actual  conditions,  such  as  finite  tunnel 
length,  transmission  of  energy  through  the  walls,  nonlinearities  at 
higher  amplitudes,  and  turbulence  in  the  flow  that  gives  rise  to  damping, 
pure  resonance  is  unobtainable.  However,  it  may  be  seen  by  examining 
figure  4(d)  that  the  lift  and  moment  were  reduced  to  20  percent  of  the 
values  away  from  resonance. 


Remarks  on  Some  Differences  Between  Theory  and  Experiment 

In  this  section,  some  limitations  and  possible  reasons  for  differ¬ 
ences  between  theory  and  experiment  are  listed  and  discussed. 

In  the  comparison  between  theory  and  experiment  shown  in  figure  4, 
an  almost  constant  difference  of  10  to  15  percent  between  the  magnitudes 
is  to  be  noted,  whereas  the  phase  angles  are  in  good  agreement.  Of  the 
several  possible  reasons  for  these  differences  between  theory  and  exper¬ 
iment,  perhaps  the  more  important  ones  axe  airfoil  thickness,  Reynolds 
number,  finite  tunnel  length,  transmission  of  energy  through  walls, 
dissipation  of  the  pressure  waves  due  to  turbulence,  and  integration 
procedures.  The  effect  of  all  these  possibilities  is  not  known  for  the 
oscillatory  case.  In  the  steady-state  case,  however,  it  is  known  that 
the  effect  of  increasing  thickness  is  to  increase  the  slope  of  the  lift 
curve. 

In  the  considerations  of  the  analytical  integration,  a  collocation 
scheme  was  used  to  solve  the  integral  equation.  In  general,  the  three 
collocation  points  used  were  found  to  be  satisfactory  as  pointed  out 
previously. 

Twelve  pressure  cells  were  used  for  the  experimental  integration  of 
the  forces.  A  stepwise  integration  procedure  was  employed;  that  is,  the 
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pressure  as  recorded  by  the  pressure  cell  was  multiplied  by  an  area  of 
the  wing  over  which  it  is  assumed  that  the  cell  will  give  an  average 
pressure.  This  procedure  gives  good  results  except  possibly  at  the 
leading  edge  where  the  pressure  variation  is  very  great.  Theoretically 

the  pressure  approaches  infinity  at  the  leading  edge  (as  ~  and 

\  \fx/ 

experimentally  it  is  found  to  be  very  large.  As  a  matter  of  fact,  if  a 
theoretical  distribution  of  pressure  is  assumed  to  be  cot  &  = 
as  shown  in  the  following  sketch  and  this  curve  is  integrated  in  the 


x 


/ 


1  +  x 


Illustration  of  the  differences  in  integration 


same  manner  as  the  experimental  curve  was  integrated  -  that  is,  by  calcu¬ 
lating  the  ordinate  at  the  same  value  of  the  abscissa  at  which  the 
pressure  cells  were  located  for  the  experiment  -  it  is  found  that  the 
area  as  determined  by  the  approximate  method  is  8  percent  less  than  the 


area  as  determined  by  integrating 


in  closed  form. 


It  is  apparent 


that  the  neglected  area  (shaded)  can  be  appreciable.  In  the  actual  exper¬ 
iment,  in  which  a  fairly  thick  airfoil  was  used,  the  neglected  area  would 
probably  be  smaller  but  could  perhaps  contribute  to  the  almost  constant 
difference  in  magnitudes  of  lift  and  moment  between  the  theoretical  and 
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experimental  results.  The  effect  of  this  integration  difference  on  the 
phase  angles,  which  have  been  shown  to  be  in  good  agreement,  would  not 
be  as  pronounced. 


Effect  of  a  Variation  in  Mach  Number  at  Constant  Tunnel  Height 


An  analytical  investigation  has  been  made  of  the  effects  of  a  vari¬ 
ation  in  Mach  number  at  constant  tunnel  height  on  the  forces  on  an 
oscillating  wing.  Some  of  the  results  of  the  previous  section  are 
employed  together  with  results  of  additional  calculations.  The  magni¬ 
tudes  and  associated  phase  angles  of  both  the  lift  force  and  the  moment 
have  been  determined  for  a  pitching  wing  and  also  for  a  wing  undergoing 
vertical  translation.  Calculations  have  been  made  for  a  constant  value 
of  the  height- semi chord  ratio  H  of  7*60  and  for  Mach  numbers  of  0.3* 
0.5>  0.7>  and  0.8.  Results  of  the  calculations  are  shown  in  figures  5 
and  6  for  the  lift  and  in  figures  7  and  8  for  the  moment. 


The  magnitudes  of  the  forces  and  moments  are  presented,  as  in  the 
previous  section,  in  the  form  of  ratios:  |La/LaT|,  J M^/M^ 1 J , 
and  l^h/^h'l •  The  phase  angles  related  to  these  ratios  are  presented 


as  a  difference  between  a  wing  in  a  tunnel  and  a  wing  in  free  air.  The 
magnitudes  and  phase  angles  are  plotted  against  a  frequency  parameter 
o)H/a,  where  o>  is  the  circular  frequency  of  oscillation  of  the  wing, 

H  is  the  height  of  the  tunnel,  and  a  is  the  velocity  of  sound.  At  a 
particular  value  of  the  frequency  parameter,  a  progressively  larger  effect 
of  the  walls  is  indicated  as  the  Mach  number  increases .  At  all  Mach  num¬ 
bers,  the  lift  is  reduced  to  zero  at  the  resonant  frequency.  The  dip  in 
the  curves  against  frequency  ratio,  which  appears  to  be  characteristic 
of  the  low  Mach  number  cases,  gradually  disappears  as  the  Mach  number  is 
increased.  As  in  the  case  of  the  magnitude  of  the  lift,  the  effect  of 
the  walls  on  the  phase  angle  increases  as  the  resonant  frequency  of  the 
tunnel  is  approached  and  also  as  the  Mach  number  is  increased. 


The  magnitude  of  the  moment  about  the  midchord  is  shown  in  figures  7 
and  8.  The  curves  have  the  same  shape  as  the  lift-ratio  curves  and  again 
decrease  to  zero  at  the  resonant  condition.  The  corresponding  phase 
angles  are  shown  in  these  figures.  Note  that  in  figures  5  to  8  only 
slight  differences  appear  between  results  for  pitch  and  those  for 
translation. 


Effect  of  Variation  in  Ratio  of  Tunnel  Height  to  Wing  Semichord 

The  effect  on  the  lift-force  ratio  and  on  the  associated  phase  angles 
of  varying  the  ratio  of  tunnel  height  to  wing  semichord  is  illustrated  in 
figure  9.  The  results  presented  in  figures  7  8  have  been  based  on  the 
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consideration  of  a  distribution  of  pressure  doublets  over  the  chord  of 
the  airfoil  and  on  satisfying  the  downwash  condition  at  three  chordwlse 
stations.  Results  for  figure  9  have  been  obtained  by  the  simplified 
procedure  of  concentrating  the  loading  in.  a  single  doublet  at  the 
quarter  chord  and  satisfying  the  downwash  condition  at  the  three-quarter 
chord.  This  approach  gives  fairly  good  agreement  with  the  results  of 
the  more  elaborate  procedure  except  near  the  resonant  frequency. 

Calculations  have  again  been  made  for  an  airfoil  oscillating  in 
pitch  about  its  midchord  for  values  of  the  height- semi chord  ratio  H 
of  7.60,  16,  and  32  at  M  =  0.3  and  0.8.  In  figure  9(a)  the  lift 
ratio  (Ia/Io/I  is  plotted  against  the  reduced-frequency  param¬ 
eter  k  =  bco/U.  Plots  for  both  Mach  numbers  are  made  to  the  same  scale 
for  ease  of  comparison.  It  is  again  apparent  that  the  effect  of 
reducing  the  Mach  number  is  to  reduce  the  effect  of  the  tunnel  vails 
and  to  raise  the  value  of  the  critical  frequency  at  which  resonance 
can  occur  for  a  given  tunnel.  For  example,  for  H  =  7*^0,  at  M  =  0.8, 
the  critical  value  of  k  is  0.30,  whereas  for  M  =  0.3*  the  critical 
value  of  k  is  increased  to  1.31*  Also,  as  was  to  be  expected, 
increasing  the  height  of  the  tunnel  has  a  marked  effect  in  reducing 
the  influence  of  the  tunnel  walls  for  most  of  the  frequency  range. 
However,  the  critical  frequency  is  reduced  for  increased  tunnel  height 
so  that  in  the  large  tunnel  the  range  of  k  below  the  fundamental 
resonance  becomes  smaller.  This  reduction  in  frequency  would  seem  to 
be  a  disadvantage  of  the  larger  tunnels.  However,  the  second  branch  of 
the  curve  for  H  =  32  at  M  =  0.3  shows  that,  for  frequencies  between 
the  first  and  second  resonant  points,  the  effect  of  the  walls  on  the 
magnitude  of  the  lift  is  no  greater  than  for  frequencies  below  first 
resonance.  Note  also  that  the  approach  to  resonance  is  quite  abrupt; 
consequently,  only  a  small  range  of  frequencies  very  close  to  resonance 
would  be  critical  and  tests  could  then  be  conducted  between  the  critical 
frequencies . 

In  figure  9(^0  the  phase  angles,  associated  with  the  results  of 
figure  9(a),  for  both  the  wing  in  a  tunnel  and  for  a  wing  in  free  air 
are  shown  as  a  function  of  reduced-frequency  parameter  k  for  values 
of  height- semi chord  ratio  H  of  7.60,  l6,  and  32  at  M  =  0.3  and  0.8. 

At  M  =  0.3,  the  effect  of  walls  on  the  phase  angle  is  generally 
very  small;  at  M  =  0.8,  the  effect  is  small  at  low  frequencies  but 
increases  greatly  as  the  critical  value  of  k  is  approached.  As  the 
resonant  frequency  is  approached,  results  at  both  Mach  numbers  show  that 
the  phase  angle  increases  in  negative  value  and  appears  to  approach  -90°. 
As  the  resonant  condition  is  exceeded,  there  is  a  sudden  shift  in  phase 
angle.  This  change  is  similar  to  that  found  for  the  oscillation  of  a 
simple  undamped  spring -mass  system  where  an  abrupt  change  in  phase  angle 
of  l8o°  is  found  as  the  resonant  frequency  is  exceeded.  Because  of  the 
complexity  of  the  kernel,  which  involves  an  infinite  series  of  Hankel 
functions,  the  phase  angle  at  resonance  has  not  been  determined. 
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In  figure  10,  the  plot  from  figure  9(a)  for  M  =  0-5  and.  H  =  7-60 
is  compared  with  some  results  of  Reissner  (ref.  2)  for  the  effect  of 
walls  on  the  lift-force  ratio  in  incompressible  flow.  This  curve  for 
M  =  0  significantly  duplicates  the  rather  large  wall  effect  at  low 
values  of  k  which  is  noted  for  the  M  =  0.3  result.  (At  low  values 
of  k,  the  curves  for  the  two  different  Mach  numbers  are  almost  coinci¬ 
dent  probably  because  of  the  fact  that  the  slightly  lower  height- 
semichord  ratio  H  =  6.28  used  by  Reissner  counteracts  the  effect  of 
decrease  in  Mach  number.)  For  values  of  k  greater  than  0.5,  the 
curves  for  the  two  Mach  numbers  separate)  the  lift-force  ratio  at 
M  =  0  approaches  unity  and  exhibits  no  effects  of  resonance  because 
the  resonant  phenomenon  arises  only  from  the  effects  of  compressibility. 


CONCLUDING  REMARKS 


This  paper  has  dealt  with  the  problem  of  an  airfoil  oscillating 
between  plane  walls  in  subsonic  compressible  flow.  It  constitutes  a 
continuation  of  the  work  initiated  in  MCA  Rep.  1150  in  that  a  method 
of  solving  the  integral  equation  is  presented  and  some  experimental 
results  are  compared  with  theory. 

The  comparison  between  theory  and  experiment  for  the  phase  angles 
between  lift  force  or  moment  and  position  is  shown  to  be  very  good, 
whereas  the  comparison  between  theory  and  experiment  for  the  magnitudes 
of  the  lift  and  moment  is  not  as  good)  however,  the  trends  are  all 
accurately  predicted.  In  all  cases  the  resonant  frequency  was 
accurately  predicted.  The  cause  of  the  apparent  discrepancy  between 
the  theoretical  and  experimental  lift  and  moment  may  be  attributed  to 
several  factors  such  as  dissipation  of  the  pressure  waves  due  to  turbu¬ 
lence  of  the  air  flow  and  transmission  of  the  energy  through  the  tunnel 
walls.  In  addition,  the  theoretical  work  was  based  on  the  concept  of 
a  very  thin  wing  at  infinite  Reynolds  number ,  whereas  the  experiments 
were  made  with  a  10-percent-thick  wing  at  a  Reynolds  number  of  approxi¬ 
mately  5  X  lO^.  The  effect  of  thickness  and  Reynolds  number  have  not, 
as  yet,  been  delineated  for  the  oscillating  case. 

As  would  be  expected,  it  is  shown  theoretically  that  the  larger 
the  tunnel  the  less  the  effect  of  the  walls.  The  critical  frequency, 
however,  is  also  reduced  as  the  tunnel  height  is  increased,  but  it  is 
shown  that  tests  may  be  made  above  the  resonant  frequency  with  no  larger 
tunnel-wall  effect  than  is  found  below  the  resonance.  In  addition,  the 
range  of  influence  of  the  resonant  region  is  greatly  reduced  so  that  only 
a  small  range  of  frequencies  need  be  avoided.  Wall  effects  are  shown 
theoretically  to  be  more  pronounced  as  the  Mach  number  is  increased  and 
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at  high  Mach  numbers  are  found  to  be  large  even  at  frequencies  well 
removed  from  resonance. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  January  12,  1955* 
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APPENDIX 

REDUCTION  OF  INTEGRAL  EQUATION  TO  THE  CASE  OF  ZERO  FREQUENCY 


In  this  appendix,  the  integral  equation  for  the  downwash  for  a  wing 
oscillating  in  a  compressible  medium  in  the  presence  of  wind-tunnel  walls 
is  reduced  to  the  zero-frequency  condition. 


If  equation  (l)  of  the  text  is  written  as 


w(x) 


lim  -&R  f  L(x0)  [aJC(M,  z)  +  aK(M,z,H)~|  dx0 

03 — >0  PU  J-1  L  J 


(Al) 


and  the  limit  taken  as  03— >0,  it  will  he  found  that  all  the  terms  of 
oiC(M,z)  and  cdK(M, z,H)  vanish  except  terms  involving  Hi(2) .  These 
terms  become  infinite;  however,  as  cd — >0,  the  asymptotic  expansion 
for  very  small  values  of  the  argument  may  be  used.  Therefore, 


Hi(2)  (nRn) 


2 

jtijjtRn 


and 


iz  P 

r2  (o\,  M(x  “  xo)  -2Map  (x  -  xD) 

lim  cseP  Hi' 2)  (p.Rn) i  - - -  =  — - - - — 


o>-*0  jt[(x  -  x0)2  +  p2(nH)2J 

The  vertical  induced  velocity  may  then  be  written  as 


w(x)  '  ^ 


h,!l  (-1)” 


X  -  X, 


X  -  X 


O  n=l 


(x  -  x0)2  +  p2(nH)2 


dx. 


(A2) 


or 


w(x)  =- 


Mab  f1  . 

h L(Xo) 


+  2 


— (X  -  x0)  n=l 

PH 


(-Dn 


rt(x  -  x0) 
EH 


Jt2(x  -  X0)2  2  2 

— ^ - 21—  +  n^Jt 


p2H2 


dx. 
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Equation  (A3)  may  be  written  as 


w(x)  =  Ztob- 
2pU^H 


f  L(xo) 

[csch  ”(X  -  *°>] 

'-1 

PH 

dx 


(A4) 


The  additional  induced  velocity  due  to  the  presence  of  tunnel  walls 
for  the  steady-state  case  in  compressible  flow  is  given  by  equation  (40) 
of  reference  13.  Equation  (A2)  can  be  reduced  to  the  same  form  by  making 
the  approximation  that  the  airfoil  chord  is  small  compared  with  the  tunnel 
height. 
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Figure  2.~  Schematic  drawing  of  model  and  oscillating  mechanism. 
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Figure  3*-  Schematic  diagram  of  instrumentation. 
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Figure  A.-  Continued. 
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(d)  M  =  0.7. 
Figure  4.-  Concluded 


Figure  5-~  Effect  of  variation  in  Mach  number  at  constant  height- 
semi  chord  ratio  H  =  7.60  on  the  magnitude  and  phase  angle  of 
the  lift  of  a  pitching  wing. 


Figure  6.-  Effect  of  variation  in  Mach  number  at  constant  height- 
semi  chord  ratio  H  =  7.60  on  the  magnitude  and  phase  angle  of 
the  lift  of  a  translating  wing. 


Figure  7*-  Effect  of  variation  in  Mach  number  at  constant  height- 
semi  chord  ratio  H  =  7*60  on  the  magnitude  and  phase  angle  of 
the  moment  of  a  pitching  wing. 
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Figure  8.-  Effect  of  variation  in  Mach  number  at  constant  height- 
semichord  ratio  H  =  7-60  on  the  magnitude  and  phase  angle  of 
the  moment  of  a  translating  wing. 
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Figure  9»-  Effect'  of  a  variation  in  tunnel  height  on  the  lift  of 
pitching  wing  for  M  =  0.5  and  M  =  0.8. 


Figure  9.-  Concluded. 


Figure  10.-  Comparison  of.  lift  ratio  for 


CM  1-1  QJ  csi 

TT  05  g<  05*  .  ”  2 

•g  >>  c  a)  CD 

“  £  £  O  o  3 

gHg*°  .S 
§  o  3  *  e  H  g 
H  H  o  c  3  C  ^ 
,  T3  g  2  £.2  «  ◄ 
c  a  h  c  °  5  u 

£  0)  O  §  O  *3  < 

^KOOSSSz 
«  «  «  a  S  ci 


fi  (D  CD 
i2  .co 

«"*  g 

dT  -2  ^ 

rt  ©  <u  •<! 

5*  o  «  U 

£  I  § 


H-asa 


i§2 

w  ^  m  J 

o  w  z*  a 

a  >  •i  p 

§S3£ 

2<  j  < 

ss»g 

s|Ib 

©  a  h  3 

IS&8 

Q  Q  W  5 

&<  ft  fc 

»  oj  wO 

H  ®<JUC(1 

g^s^« 

H  g|og 

oS  w  ft  w 

<  rt  K  S  W 

S5ZHHH 


T3  O 

«§I 

d  «  a 

S3  r 

w  «  in 

sli 

p«  **  5 
5  m  . 
o  M  o- 

u:i5 

ag  . 

2  Q  IO 

ft  H  irt 


S  R§ 

J  ft;  ^ 

3  * 

O  ^  “S' 

PJz 

•g  a  ts  z 

?sl< 

|3“s 


,  rt  ® 

1  o  *3 

I  ©  T-*  S*  ««-( 

1  %  a~  ° 

:  a  a  ®  a 
•2<stJ 

■a  .  ■"  .2 
rt  §  °  ® 
T3  ©  «  X 

©  X  ®  73 

rti2 

©  P  g  «♦-. 
t*  o  3  o 

J  ®  o  ®  ’S 

->  P 
©  *»  °  ^ 
i'g  ©  § 

3  X  .q 

;§•§  s 
3*8  cJ 

1  M  8  -  c 
.5  g  2  O . 
S  <=  ®  ® 
ta  rt  ■g  e 

1  e  c  2  5 
3  .5  S  s 


°  ■a 

ass 

rt  C  a 
l  2  © 

O  a 

*8  5  * 

rt  rt  <2 

M  g  J2 

©  o  C 

f°S 

3  .  © 
c  T3  3 

•g  S  § 

s  a  s  • 

&S32 


3c© 
£  *>  £  -a 

3.Sg§ 

o’  >  .5  3 

©  R  ^  G 

^  o  ©  - 

C  ~  &  P 

.3  .g  © 

Krv  ^ 


-pH  -*->  —  ,rM 

©eg© 
ra  ©  o  S 
o  w  12  is 


X!  -  *  S 

rt  ©  ©  O 


CM  rH  q,  CM 
•  •  3  •  03 

asfsj, 

5  >>1 

•2  S  t  | 
1® 

3  -S  a  *  g 

E-t  u  o  e  5 
rt  ©  at  oo 
'O  ©  5h  >,  72 

g  m  u  q  o 
r  ©  o  3  o 
?  BS  U  «  fe 


«aSS 


<J  TJ  O 

gzgw  S| 
,HORd  C  a 

ra  co  r.  j  ffl  o 

U  B  W  M  H  -W  r 

i23«§5  sp 

<el§8”^ 

si§^o1'r 

§  S  ^  ^  §  o  S 

IftOOojS 

laesi|i 

■sgsM- 

H  ®  <  H  w  W 

3!  >UBW^  ^  ^ 

uaw§wgoou 


«  0) 

3  £al- 

a  *e  rl  o 
;  a  a  §  n 

2<2t5 

*  T3  <2 

C  m  O  ® 

5  a  ~  © 

>  -c  ci  ®  i 

-  a  a  ®  a 
^■SiS  as 

mi's 

P.5  S  S-g 

i  ©  -*-»  °  -pj 

!  §Sjs  “ 

if  g-l  s 

^  T3  2  © 

;  e  rt  c  jt 

e  u  o  ©  i 
•HQ 

jt  ««h  q>  © 

he  rt  'g  g 

>.5  c  g  f 

2  ^  .2  *3  : 

in  i2  _  —i  . 
u  q  m 
»  3  n  © 


°  a 

.2.2  S 

^  rt  c  ?T 

l^oj 

•  a  5s 
I  «.■  fs 

®  o  «s 

l°a 

i-dg 

r-  ?h  g 

©  o  2 

j  h  ci 

So® 

s  a  a 

1  kT  g  --3  CO 
&  M  -2  ^ 

g“S^ 

§  .S  g  g 

cr  ^  .5  3 

©  ^  Jh  cs 
Jh  o  © 

H-»  £V  X3 

sals 

»§^S 

fi  ©s  O 
P  J3  >  .3 
XI  _  ?  C 

rt  ©  ©  O 
•H  C  T3  ® 
^  H  Ji  JQ 


SzzM 
C  £  o  j 


Ps 

2  H  u 


•4-J  txj 

•d  Ph 

a  x 

a  w 

o  Q 

&< 

2.5  < 

2  <  H 

t-|  S 

S  5  w 


Hh4£ 

8^1 

e*l 

li,  <  <; 

ft  Z  5 

H00 

III 

Si? 

Sag 

H  S  ? 

H  H  h 


>,  fi  % 

s'as 

M  T3  ^ 

73  H 

>  S  < 

j  .3 

fe<6 


,  rt  3 
q®5? 
rt  ^ 'd  13 
-H  aa2 

°2  <qd 

a;  a 

5  S  co  o 

Otj  o  ID 

ft  ©  r3  © 
2  ©  S3  g 

5-a  S 

c  .2  8  8 


c  s 

.2  § 

d  fi  c 


d  ?  °  ffl 

Saul'S 

bfl  C3  c  2 

©  -3  -5  S 
e  j3  t3  © 

— '  Sd  CD 

®  o  C  0 

h  §  S  a 


<N  H  ,JJ  CM 
^  05  q«  05  .  ®  2 

'■“'•a  J2  d 

X  >.g  geo 

“  ©  £  o  CD  Xj 

gH  “  .co 

§  JC  o  w  -<C  s- 

3  ©  xl  ^  S  5 
H  »h  o  q  2  kT  ^ 

_  rt  ©  rt  to  ^  ^ 

T3  ©  ti  >1-^  S  /! 

c  w  e  c  o  ,g  U 

^  ©  o  3  O  rt  < 

^«U«^S5z; 

CM  CO  hhhh 


2.§< 
**  >  n 
CO  TJ  W 
;z<  H 

Sl8 

UiW 
55  ^  Eh 


iSSS 

!MB- 

'i«s 

wgol 

lly' 

e^s( 

ftdg1 
o  o  S 

H§S 

goM 

[*<  <  < 

“OS' 

b  S  S ' 

o  ga  „ 

±t 


j  g* 

«  T3  Z 

r, 


S8|« 

>h  Qftjj 
Qr]  -*h 

3|  ©  o 
xi  h  e  *• 
o  ^  o  w 
fn  o  x  © 
a-S^  a 

®  2  -g  a 

£  il  «  § 


fj  O  Q 
.2  tD  ■*->  u 


.d  « is  a 

ti  e  rt  *3 

1*8  ^ 
®  .3  g  2 
xl  ^  © 

2  Wirt'S 
ho  q  a  g 
©  Tj  *5  S 
G  rt  XJ  w 
d3  © 

®  o  e  c 
£  ra  ©  .2 
H  o  to  X 


the  influence  on  the  magnitude  of  wall  effects  of  tiie  influence  on  the  magnitude  of  wall  effects  of 

variations  in  frequency,  Mach  number,  and  ratio  of  variations  in  frequency,  Mach  number,  and  ratio  of 

tunnel  height  to  wing  semichord.  Comparison  is  tunnel  height  to  wing  semichord.  Comparison  is 

made  with  experimental  measurements  for  several  made  with  experimental  measurements  for  several 

subsonic  Mach  numbers.  subsonic  Mach  numbers. 


